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Claims 

1. A method for identifying pathology in a brain image comprising the 
steps of: 

(a) determining the location of the midsagittal plane (MSP) by 
5 calculating up to 16 approximated fissure line segments (AFLSs) and 

removing outlier AFLSs, the outlier AFLSs having a larger angular 
deviation from the MSP than a predefined threshold; 

(b) calculating the ratio of the number of outlier AFLSs to the 
number of inlier AFLSs, the inlier AFLSs having a smaller angular 

10 deviation from the MSP than the predefined threshold; and 

(c) comparing the ratio with a further predetermined threshold 
value, the ratio exceeding the further predetermined threshold value 
when pathology is present in the brain image. 

2. A method as claimed in claim 1 , wherein the image comprises image 
15 data and the method further comprises the step of: 

reformatting the image data to align the image with the midsagittal 

plane. 

3. A method as claimed in claim 2, wherein the reformatting is performed 
using vector manipulation, 

20 4. A method as claimed in claim 2 or 3, wherein the midsagittal plane has 
an equation x'+d'=0 where x' is a co-ordinate in the OX direction and d' is a 
constant, the image having an original co-ordinate system OXYZ before 
reformatting and a reformatted co-ordinate system after reformatting 
comprising co-ordinates OX, O'Y' and OZ, the step of reformatting 

25 comprising: 
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(i) starting with an MSP equation : 
ax + by + cz +d = 0, 

where d is less than 0 and (a,b,c) is the unit normal vector of the MSP; 

(ii) determining two points that are the intersections between the MSP and 
5 the volume of the brain shown in the image under examination, denoting the 

two points as A (xa, Ya, 0) and (xe, Yb, 0), both of the two points being on the 
MSP; 

(iii) when a is not equal to 0, calculating A and B as follows: 
x A =-d/a, y A =0 

10 x B = -(d+b(ySize -1 ))/a, y B =ySize - 1 

where ySize is the number of voxels in the Y direction of the original volume; 

(iv) when a is 0, b should not be 0, and calculating A and B as follows: 
x A =0, y A = -d/b 

xb = x Size -1 , y B = -(d + a(xSize -1 ))/b 

15 where xSize is the number of voxels in the X direction of the original volume; 
and 

(v) changing the original coordinate system OXYZ to the reformatted co- 
ordinate system O' X 1 Y'Z, the unit vectors of the co-ordinates 0'X\ O'Y' and 
O'Z' being calculated in the following manner: 

20 O' X' = (a, b, c) = (n x i, n y1t n z1 ) 

CTY' = ((xA-xB)/ A-B , (yA - yB)/ A-B ,0) = (n^, ^ n^) 
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O'Z' = O'X'xO'Y' = (n x3 , n y3 , n z3 ) 

where A-B is the Euclidean distance between points A and B, the 
transformation between OXYZ and O'X'Y'Z' being defined as follows: 

X' = n x i X + n y1 Y + n z i Z 

5 Y' = n x2 X + n y2 Y + n^Z 

Z' = n x3 X + n y3 Y + n z3 Z 

O' = 0. 

5. A method according to any one of the preceding claims further 
comprising the step of: 

10 processing one or more slices through the image to determine 

interhemispheric symmetry. 

6. A method according to claim 5, wherein the step of processing 
comprises processing one or more slices having a plane substantially 
perpendicular to the midsagittal plane. 

15 7. A method according to claim 5 or 6 wherein the step of processing 
comprises processing one more slices having a plane which intersects the 
midsagittal plane. 

8. A method according to any one of claims 5 to 7 wherein the image 
comprises data representing the brain and background image data, the step 

20 of processing comprising extracting the background image data and retaining 
the data representing the brain. 

9. A method according to claim 8, wherein the image is comprised of a 
number of pixels, the step of extracting comprising calculating a histogram of 
grey level of the image along a first axis against the number of pixels of the 
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image along a second axis, locating a first peak in the histogram formed at the 
intersection of a positive slope and a negative slope, extrapolating the 
approximating the positive slope back to intersect the first axis and taking the 
value of grey level at the intersection of the positive slope with the first axis as 
5 a threshold grey level value. 

10. A method according to claim 9, wherein the step of processing 
comprises further processing the image to obtain a binary image by setting 
the pixels having a grey level greater than the threshold grey level to 1 and 
the pixels having a grey level lower than the threshold grey value to 0. 

10 11. A method according to claim 10, wherein the image depicts the brain 
and a skull, the method further comprising the steps of opening the binary 
image with a structuring element of pixels and expanding the structuring 
element to fill the image of the skull to produce a head mask. 

12. A method according to claim 1 1 , further comprising cropping the image 
15 to the head mask to remove the background data visible in the image. 

13. A method according to claim 12, further comprising calculating a 
histogram of grey level of the image along a first axis against the number of 
pixels of the image along a second axis for each hemisphere, smoothing the 
histogram for each hemisphere, comparing the histograms of the two 

20 hemispheres by calculating the difference between the histograms. 

14. A method according to claim 13, wherein the step of calculating the 
difference in values between the histograms comprises calculating a 
difference function between the values of the histograms for all grey levels. 

15. A method according to claim 13, wherein the difference function has a 
25 maximum value and a minimum value, the method further comprising 

identifying when the difference between the maximum value and the minimum 
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value exceeds a threshold difference value to indicate pathology presence in 
the brain under examination. 

16. A method according to any one of claims 13 to 15, wherein the step of 
comparing the histograms comprises calculating the difference in area of the 

5 histograms, or the normalized difference in area of the histograms. 

17. A method according to any one of claims 13 to 16, further comprising 
the step of aligning the histograms along an axis before comparing the 
histograms. 

18. A method according to claim 17, wherein the step of aligning the 
10 histograms comprises maximising the cross-correlation of the histograms. 

19. A method according to any one of claims 6 to 18, further comprising 
determining extent of pathology in the brain under examination. 

20. A method according to claim 19 wherein the step of determining extent 
of pathology comprises determining the ratio of the number of slices with 

1 5 pathology to the total number of slices studied. 

21. A method according to claim 19 or 20, further comprising determining 
the location of the pathology. 

22. A method according to claim 21, wherein the step of determining the 
location of the pathology comprises locating the hemisphere that contains 

20 pathology. 

23. A method according to claim 22, wherein the step of locating comprises 
determining the modality of histograms for both hemispheres of the brain 
under consideration, tri-modality existing for both hemispheres being an 
indication of a small pathology existing in one or both hemispheres. 
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24. A method according to claim 23, comprising determining if one or both 
of the hemispheres is not tri-modal indicating pathology existence therein and, 
if one or both of the hemispheres is not tri-modal determining the modality of 
the histograms for each slice. 

5 25. A method according to claim 24, further determining if tri-modality 
exists for all of the slices within both hemispheres, and if it does not exist 
determining if the ratio of the number of outlier AFLSs to the number of inlier 
AFLSs. 

26. A method according to claim 24, further determining if tri-modality 
10 exists for all of the slices within both hemispheres, and if it does not exist 

determining if the ratio of the number of outlier AFLSs to the number of inlier 
AFLSs. 

27. A method according to any one of the preceding claims, further 
comprising the step of calculating one or more additional features of each 

15 hemisphere, comparing the one or more additional features, and if their 
difference is larger than a predetermined threshold, signalling the existence of 
pathology. 

28. A method according to claim 27, wherein the step of calculating one or 
more features comprises calculating one or more of intensity mean, moments, 

20 cross-correlation, and mutual information. 

29. An apparatus arranged to perform a method for identifying pathology in 
a brain image according to any one of the preceding claims. 

30. A computer program product comprising computer program 
instructions readable by a computer apparatus to cause the computer 

25 apparatus to perform a method according to any one of claims 1 to 28. 



